Field dependent quasiparticles in the infinite dimensional Hubbard model 
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We present dynamical mean field theory (DMFT) results for the local spectral densities of the one- 
and two-particle response functions for the infinite dimensional Hubbard model in a magnetic field. 
We look at the different regimes corresponding to half-filling, near half-filling and well away from 
half-filling, for intermediate and strong values of the local interaction U . The low energy results are 
analyzed in terms of quasiparticles with field dependent parameters. The renormalized parameters 
are determined by two different methods, both based on numerical renormalization group (NRG) 
calculations, and we find good agreement. Away from half-filling the quasiparticle weights, Za{H), 
differ according to the spin type a =| or a =],. Using the renormalized parameters, we show that 
DMFT-NRG results for the local longitudinal and transverse dynamic spin susceptibilities in an 
arbitrary field can be understood in terms of repeated scattering of these quasiparticles. We also 
check Luttinger's theorem for the Hubbard model and find it to be satisfied in all parameter regimes 
and for all values of the magnetic field. 

PACS numbers: 72.15.Qm, 75.20Hr, 73.21.La 
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I. INTRODUCTION 

A feature of strongly correlated electron systems, such 
as heavy fermions, is their sensitivity to an applied mag- 
netic field, which makes a magnetic field a useful exper- 
imental probe of strong correlation behavior. A man- 
ifestation of this sensitivity is the very large paramag- 
netic susceptibility observed in these systems. In terms 
of Fermi liquid theory, the large paramagnetic suscep- 
tibility can be interpreted as due to quasiparticles with 
exceptionally large effective masses. These large effective 
masses arise from the scattering of the electrons with the 
enhanced spin fiuctuations induced by the strong local 
Coulomb interactions. An applied magnetic field sup- 
presses the spin fiuctuations causing a reduction in the 
effective masses, which can be seen experimentally in de 
Haas- van Alphen measurementsi"^. Not only do the effec- 
tive masses depend on the magnetic field, they may also 
differ for the spin up and spin down electrons^i^. Another 
feature that refiects the enhanced sensitivity to an ap- 
plied field is metamagnetic behavior, where the spin sus- 
ceptibility x{H) in a finite field H increases with the field 
strength such that dx{H) / dH > 0, which has been ob- 
served in some heavy fermion compounds^. This can be 
understood in terms of an increase in the effective mass 
for larger fields opposite to the effect described above. 
It is related to the fact that strong magnetic fields can 
induce localization in narrow conduction bands as pre- 
dicted theoretically^"^. This has been observed experi- 
mentally, for instance, in quasi-two dimensional organic 
conductors^.. 

A lattice model which can mimic many of these effects 
is the single band Hubbard model. It has played a similar 
role for lattice models as the Anderson model for impu- 
rity models, being the simplest model of its type, where 
the interplay of kinetic energy and strong local interac- 
tions can be studied. Here we are interested in studying 
the magnetic response of this model for different inter- 



actions and fillings. The calculations are based on the 
dynamical mean field theory (DMFT) with the numer- 
ical renormalization group (NRG) to solve the effective 
impurity problem. We present results for the local spec- 
tral densities and spin dynamics in parameter regimes 
at half filling and for finite hole doping, where different 
responses to the magnetic field can be observed. 

We interpret the low energy results in terms of quasi- 
particles which are characterized by field dependent 
renormalized parameters. The approach is similar to 
that used earlier-iSiii for the Anderson model in a mag- 
netic field. We deduce the renormalized parameters by 
two different methods based on the NRG calculations. 
These parameters can be used to define a free quasipar- 
ticle density of states, which gives the asymptotically ex- 
act spectral behavior for low energy. The effects of the 
interactions between these quasiparticles can be taken 
into account using a renormalized perturbation theory 
(RPT) ^^i^^ . It was shown earlier that a very good de- 
scription of the r = spin and charge dynamics for the 
Anderson model in the Kondo regime can be obtained 
by summing the RPT diagrams for repeated quasipar- 
ticle scatteringii. Here we extend these RPT calcula- 
tions for the spin dynamics to the lattice case and show 
that we can also understand the DMFT-NRG for the lo- 
cal dynamic spin susceptibilities in terms of quasiparticle 
scattering. 



II. DYNAMICAL MEAN FIELD APPROACH 
AND RENORMALIZED PARAMETERS 

The Hamiltonian for the Hubbard model in a magnetic 
field in the grand canonical formulation is given by 
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where tij are the hopping matrix elements between sites 
i and j, U is the on-site interaction; /io- = /i + ah, where 
/U is the chemical potential of the interacting system, and 
the Zeeman splitting term with external magnetic field H 
is given by h — g^BH/2, where /j,b is the Bohr magneton. 
We are dealing with the one s-band Hubbard model here, 
so no coupling of the field to angular momentum states 
has to be included. 

From Dyson's equation, the Green's function Gk,a-i^^) 
can be expressed in the form. 
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where T,a{k,uj) is the proper self-energy, and e{k) = 
e~^'''^^~^^Hij . The simplification that occurs for the 
model in the infinite dimensional limit is that Ttcy{k,uj) 
becomes a function of u onl y^^i^^ . In this case the local 
Green's function G^°'^{uj) can be expressed in the form. 
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where D{e) is the density of states for the non-interacting 
model {U — 0). It is possible to convert this lattice prob- 
lem into an effective impurity problemi^. We introduce 
the dynamical Weiss field Q^^yii-^) and write the Green's 
function in the form 
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which is equivalent to 



(4) 
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The Green's function G^°''{uj) can be identified with the 
Green's function Ga{ijj) of an effective Anderson model, 
by re-expressing Q^]^ (uj) as 



so that 



UJ — 
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with edcr = — Mct- The function Kfj{uj) plays the role of 
a dynamical mean field describing the effective medium 
surrounding the impurity. In the impurity case in the 
wide band limit we have Kc[ijj) — —iA. Here, as can 
be seen from equations ([5]) and ([4]), K^^iuj) is a function 
of the self-energy I]cr(w), and hence depends on a. As 
this self-energy is identified with the impurity self-energy, 
which will depend on the form taken for Ka{uj), it is clear 
that this quantity has to be calculated self-consistently. 
Starting from an initial form for Ka-{uj), is calcu- 

lated using an appropriate 'impurity solver' from which 
G^°^{uj) can be calculated using equation ([3]), and a new 
result for Ka-{uj) from equations ^ and ((6]). This Kduj) 



serves as an input for the effective impurity problem and 
the process is repeated until it converges to give a self- 
consistent solution. These equations constitute the dy- 
namical mean field theory (DMFT), and further details 
can be found in the review article of Georges et a\M-. 

We need to specify the density of states -D(w) of the 
non-interacting infinite dimensional model, which is usu- 
ally taken to be either for a tight-binding hypercubic or 
Bethe lattice. Here we take the semi-elliptical form cor- 
responding to a Bethe lattice, 

Die) = -^VD^^¥TiIor (8) 

where 2D is the band width, with D = 2t for the Hub- 
bard model, and fiQ the chemical potential of the free 
electrons. We choose this form, rather than the Gaus- 
sian density of states of the hypercubic lattice, as it has 
a finite bandwidth. 

Before considering in detail the methods of solving 
these equations, we look at the form of these equations 
in the low energy regime, where we can give them an in- 
terpretation in terms of renormalized quasiparticles. We 
assume that we can expand Scr(w) in powers of to for 
small Lu, and retain terms to first order in oj only. Sub- 
stituting this expansion into the equation for the local 
Green's function gives 
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and = 1/[1-E;(0)]. (10) 



G'o°:A^) = z,Jde 

where 

Ao.cr = z^ip^~j:^{o)), 

We have assumed the Luttinger result that the imaginary 
part of the self-energy vanishes at uj — 0. As the Green's 
function in equation ([9]) has the same form as in the 
non-interacting system, apart from the weight factor Za-, 
we can use it to define a free quasiparticle propagator. 
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We then interpret Za as the quasiparticle weight. We 
will refer to the density of states /5o,(j(w) derived from 
this Green's function via /5o,ct(w) = — ImGo,<T(w -t- i5)/7r 
as the free quasiparticle density of states (DOS). For the 
Bethe lattice, the quasiparticle DOS takes the form of a 
band with renormalized parameters. 



Dl - (cJ + Ao.a) 
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where Da = z^D. We can also define a quasiparticle oc- 
cupation number ri° by integrating this density of states 
up to the Fermi level. 
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With a generalization of Luttinger's theorenii^ for each 
spin component it is possible to relate this free quasi- 
particle occupation number n° to the expectation value 
of the occupation number n^y in the interacting system 
at T = 0. Using the quasiparticle density of states in 
equation lfT2|) . we can rewrite equation (flS]) as 



OO 
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where 0{e) is the Heaviside step function and D{e) as 
given in equation Assuming Luttinger's result for 
each spin component, the right-hand side of equation lfT4|) 
is equal to n„. We then have the result, 



(15) 



that the occupation for electrons of spin a is equal to the 
number of free quasiparticle of spin cr, as calculated from 
equation (flS]) . It should be noted that there is no simple 
generalization of the h = result, /j- — /j,o = E(0), in the 
spin polarized case. 

To solve the DMFT equations we need an 'impurity 
solver', and the most commonly used are the quantum 
Monte Carlo, the exact diagonalization (ED) method and 
the NRG, all of which have advantages and disadvan- 
tages. Here we use the NRG approach as it is the most 
accurate method for calculations at T = and for the 
low energy excitations. There has been a DMFT study 
of the static properties of a half-filled Hubbard model in 
a magnetic field using the ED method by Laloux et alfZ.. 
The focus of our paper here, however, is rather different 
so there is little overlap with this earlier work but, where 
there is, we make comment and compare with their re- 
sults. 

To evaluate the renormalized parameters, Za- and /io.o-, 
which specify the form of the quasiparticle DOS, we use 
two different methods based on the NRG approach. The 
first method is a direct one, where we use the NRG deter- 
mined the self-energy Scr(w) and the chemical potential 
Ha-, and then substitute into equation (fTO| for z^r and 
lio,a- The second method is indirect, and makes no ref- 
erence to the self-energy. It is based on the quasiparticle 
interpretation of the NRG low energy fixed point of the 
effective impurity. To explain this approach we need to 
consider in a little detail how the NRG calculations are 
carried out. 

In the NRG approachi^ the conduction band is log- 
arithmically discretized and the model then converted 
into the form of a one dimensional tight binding chain, 
coupled via an effective hybridization Va to the impurity 
at one end. In this representation Ka{uj) — \Va\'^go,cr{i^), 
where go^„{uj) is the one-electron Green's function for the 
first site of the isolated conduction electron chain. As 
earlier, we expand the self-energy Scr('j-') to first order in 
uj, and then substitute the result into equation ([7]). We 
then define a free quasiparticle propagator, Go.ai^), for 



the impurity site as 

Go,<t(w) = 

where 
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Ed^ = Za(ed^ + S^(0)), |KtP = 2<t|K7 
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In the DMFT approach we identify Go,(j (w) with the local 
quasiparticle Green's function for the lattice ifTTj) . 



(18) 



which specifies the form of go,<T('jj) in (flBl) and yields 
/2o,cr — ~£dcT- By fitting the low energy single parti- 
cle excitations found in the NRG results to the poles of 
this Green's function, we can deduce the parameters ida 
and Va, as has been explained in earlier work^. The 
quasiparticle weight Za- is then obtained from the re- 
lation Zcr — |Kt/VctP in equation lfT7|) . and flo,a from 

Using the DMFT-NRG approach we can calculate the 
spectral densities for the local two-particle response func- 
tions, as well as single-particle ones. The main inter- 
est here will be in the local longitudinal and trans- 
verse dynamic spin susceptibilities, x/(w) and Xti^)- 
Having calculated the renormaHzed parameters, which 
describe the free quasiparticles, we can compare the 
DMFT-NRG results for the dynamic susceptibilities with 
the corresponding quantities calculated via a renormal- 
ized RPA-like treatment that takes account of repeated 
quasiparticle-quasihole scattering. This approach has 
been described fully elsewhere for the Anderson impurity 
modeliiiii. The calculations here proceed along similar 
lines as for the effective impurity model. The equation 
for the transverse susceptibility is 



(19) 



where xui^) transverse susceptibility calculated using 
the free quasiparticles density of states given in equation 
(fT2l) . and Ut{h) is the irreducible quasiparticle interaction 
in this channel. For the Anderson model it was possible 
to calculate Ut{h) in terms of the renormalized on-site 
interaction U. Though it is possible in the lattice case to 
calculate U, we have no way of deducing Ut{h) from it 
since unlike in the impurity case we do not have an exact 
expression for Xt(0) in terms of U. We determine it to fit 
the DMFT-NRG result for Rext(a;) at the single point 
Lu = 0. The corresponding result for the longitudinal 
susceptibility xi{^) is 

= (20) 

Xtt(^,^) + Xii(^,^) +4L/;(fe)xtt(tJ,fe)xxi(t^,fe) 
2{l~AU^{h)xn{^,h)xu{^,h)) 

where the susceptibilities Xaai^) are those for the free 
quasiparticles, and Ui{h) is determined by fitting the 
DMFT-NRG result for Rex/(0). 
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Having covered the basic theory, we are now in a po- 
sition to survey the results for the Hubbard model in 
different parameter regimes. 



III. RESULTS AT HALF-FILLING 

We present the results at half-filling for three main 
parameters regimes where we find qualitatively different 
behavior. The results in all cases will be for a Bethe 
lattice with a band width W = 2D — 4, setting t — 1 
as the energy scale. The local spectral densities are cal- 
culated from equation |(3| using the NRG deduced self- 
energy^-. In the evaluation of all NRG spectra we use 
the improved metho d^^i^^ based on the complete Anders- 
Schiller basis^. In concentrating on the field induced po- 
larization, we do not include the possibility of antiferro- 
magnetic ordering. The regimes are (a) a relatively weak 
coupling regime where U is smaller than the band width, 
(b) an intermediate coupling regime with W < U < Uc, 
where Uc is the value at which a Mott-Hubbard gap de- 
velops in the absence of a magnetic field {Uc « 5.88)2^, 
and (c) a strong coupling regime with U > Uc- 



A. Weakly correlated regime 

The first plot in figure [T] gives the spectral densities for 
the majority spin electrons for various magnetic 

field values in the weakly correlated regime, U ~ 2. 

0.35 p • • • • • n 
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FIG. 1: (Color online) The local spectral density for the ma- 
jority spin Pi{lj) for (7 = 2 and various fields h. 

We can see clearly that, for increasing magnetic field, 
more and more spectral weight is shifted to lower ener- 
gies (the opposite happens for the other spin component, 
which is not displayed here). Above ~ 1.0 the system 
is completely polarized, 2m — 1. This extreme high field 
limit corresponds to a band insulator. There is a gap of 
the magnitude Ag(h) = 2h + U — W between the upper 
(minority) and lower (majority) bands, which both have 



the semi-elliptical form as for the non-interacting system 
with W = 4. At this point dynamic renormalization 
effects have vanished. The inverse of the quasiparticle 
weight Z(j(h), which corresponds to the enhancement of 
the effective mass m*^lh) — m/ z„(h)^ is shown as a func- 
tion of h in figure [21 




FIG. 2: (Color online) The inverse of the quasiparticle weight 
Za{h) calculated from the impurity fixed point (FP) and di- 
rectly from the self-energy and the magnetization m{h) for 
U = 2 and various fields h. 

We calculated the values of Za{h) using the two methods 
described earlier. At half filling we have zi{h) ~ z^{h) 
and we have plotted the average of the values for (7 = 1 
and (J — which is compared for the two methods. The 
deviation for the values for the different spins is only 
due to small numerically inaccuracies and is less than 
2%. The method based on analyzing the excitations of 
the impurity fixed point (FP) is only appHcable in the 
metaUic regime and when the system is not completely 
polarized. From these results shown in figure [2] it can 
be seen that the two sets of values are in good agree- 
ment. The values of Za{h) increase from about 0.75 to 1, 
which corresponds to a progressive "de-renormaHzation" 
of the quasiparticles with increasing field, as observed 
earlier for the impurity modeU^. Since the interaction 
term in the Hubbard model acts only for opposite spins 
it is clear that there is no renormalization when the sys- 
tem is completely polarized with one band fully occupied 
and the other empty. The expectation value of the double 
occupancy {n^ni) decreases with increasing field, which 
further illustrates why the interaction term becomes less 
important for larger fields. 

We can also follow the field dependence of the renor- 
maHzed chemical potential fio^a{h) as plotted in fig- 
ure [21 In the case of particle hole symmetry we have 
Mo.tCi) — ~Ao, 4(^)5 and similar to the case for the quasi- 
particle weight we have displayed the average of the up 
and down spin values for each method. Also here the de- 
viation within one method is only due to small numerical 
inaccuracies and less then 2%. 
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FIG. 3: (Color online) The renormalized chemical potential 
Mo,ct(^) calculated from the impurity fixed point (FP) and 
directly from the self-energy for f/ = 2 and various fields h. 



Again the agreement between the two methods of calcu- 
lation is very good over the full range of magnetic fields. 
Mean field theory is valid for very weak interactions, so 
we compare our results for /io,o-(/i) for ?7 = 2 with the 
mean field value /ijf ^ = fi + ah — t/n™^ in figure [3l The 
results coincide for h = 0, when /2™^ = and when 
the system becomes fully polarized at large field values, 
/i™^ — —a{U + h), but in general jj,™^. > fio,a{h,). We also 
compare the mean field (MF) result for the field depen- 
dence of the magnetization m{h) with the one obtained 
in the DMFT calculation in figure IH 




m(h) 

FIG. 4: The magnetization in the mean field approximation 
compared with the DMFT result for f/ = 2 and for the full 
range of magnetic fields h. 

The general behavior is similar, but the mean field theory 
without quantum fiuctuations overestimates the magne- 
tization, as one would expect. 



B. Intermediate coupling regime 

In the next plot in figure O where U — 5, we show 
typical behavior of the local spectral density in the inter- 
mediate coupling regime. 
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FIG. 5: (Color online) The local spectral density for the ma- 
jority spin p] {uj) for [/ = 5 and various fields h. 

Similar to the weak coupling regime, we find a shift of 
spectral weight towards lower energy for the majority 
spin. There is, however, a difference in the way this hap- 
pens due to the initial three peak structure, namely the 
quasiparticle peak in the middle gets narrower for in- 
creasing field and finally vanishes in the polarized phase. 
The quasiparticle weight, which is shown in figure [6l re- 
fiects this behavior by decreasing to zero with increasing 
field signaling heavy quasiparticles. Here, as in the weak 
coupling case, we plot the average of the spin up and 
down results for each method. The deviations can be 
larger here, especially close to the metamagnetic transi- 
tion. 

When the material is polarized, z^ih) reverts to 1, which 
corresponds to the band insulator as before. This ap- 
proach to the fully polarized localized state in high fields 
contrasts with that found in the weak couphng regime. 
It gives rise to metamagnetic behavior in this parameter 
regime. To ihustrate further the different response to a 
magnetic field, the real part of the local longitudinal dy- 
namic spin susceptibility x;(ti;,/i) as a function of uj is 
shown for various values of h in figure [7l 
It can be seen that the local susceptibility = 
Rexi{0,h) in this regime increases with h so that 
dx^°'^{h)/dh > 0. This can also be seen in the curvature 
of the magnetization shown in the inset of figure [H This 
is behavior characteristic of a metamagnetic transition 
and related to the magnetic field induced metal-insulator 
transition. Laloux et ali^ find metamagnetic behavior in 
a similar parameter regime. There a comparison is made 
with results from the Gutzwiller approximation, which 
gives such a behavior already for smaller values of the 
interaction, and we refer to their paper for details. 
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FIG. 6: (Color online) The inverse of the quasiparticle weight 
Zc{h) calculated from the impurity fixed point (FP) and di- 
rectly from the self-energy for U = 5 and various fields h. The 
inset shows the magnetization m{h). 
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FIG. 8: (Color online) The comparison of the spin dependent 
occupation numbers and Ucr corresponding to Luttinger's 
theorem in a magnetic field for [/ = 5 and the range of fields 
h. 
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FIG. 7: (Color online) The real part of the local longitudinal 
dynamic spin susceptibility for (7 = 5 and various fields h. 



We can also check the Luttinger theorem in a magnetic 
field ifTS]) by comparing the values of deduced from 
integrating the quasiparticle density of states, virith the 
value of Tie calculated from the direct NRG evaluation in 
the ground state. The results are shown in figure [HI 
It can be seen that there is excellent agreement between 
the results of these two different calculations, = rio-, 
so that Luttinger's theorem is satisfied for all values of 
the magnetic field in this intermediate coupling regime. 

Having deduced the renormalized parameters of the 
quasiparticles from the NRG results, we are now in a po- 
sition to test how well we can describe the low energy 
dynamics of this model in a magnetic field in terms of a 
renormalized perturbation theory. It is of interest first 
of all to see how the free quasiparticle density of states 
Pi^^aiuj) multiplied by zdK) compares with the full spec- 
tral density Pa(y>)- In figure [9] (upper panel) we make a 



comparison in the zero magnetic field case. 
We see that the quasiparticle band gives a good represen- 
tation of the low energy peak in Paiijj) and, as expected, 
does not reproduce the high energy features. These, how- 
ever, to a fair approximation can be described by the 
mean field solution pmf (w) weighted with a factor 1 — 
as can be seen in figure [9] (upper panel) . A case with 
a finite magnetic field h — 0.15, where the peaks in the 
density of states of the two spin species are shifted due 
to the induced polarization relative to the Fermi level, is 
shown in figure[9] (lower panel). The figure focuses on the 
region at the Fermi level and one can see that the free 
quasiparticle density of states describes well the form of 
Pcr{(^) in the immediate vicinity of the Fermi level. It is to 
be expected that the frequency range for this agreement 
can be extended if self-energy corrections are included in 
the quasiparticle density of states using the renormalized 
perturbation theory as shown in the impurity case^^. 

We now compare the NRG results for the longitudinal 
and transverse local dynamic spin susceptibilities for the 
same value U — 5 and a similar range of magnetic field 
values with those based on the RPT formulae ((20|) and 
(fT9|) . In figure [lO] (upper panel) we show the imaginary 
part of the transverse spin susceptibility calculated with 
the two different methods. 

It can be seen that RPT formula gives the overall form 
of the NRG results, and precisely fits the gradient of the 
NRG curve at oj — 0. Some of the relatively small dif- 
ferences between the results might be attributed to the 
broadening factor used in the NRG results which gives a 
slower fall off with to in the higher frequency range, and 
a slightly reduced peak. We get similar good agreement 
between the two sets of results for the same quantity for 
the case with a magnetic field h — 0.15, shown in figure 
[TOl (lower panel) . 

In figure [TTl where we give both the real and imagi- 
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FIG. 9: (Color online) The free quasiparticle density of states 
in comparison with interacting local spectral density for U = 
5 and h = (upper panel). We have also plotted a thin 
black line for pmf(c^) = [D{uj + U/2) + D{uj - C//2)]/2 which 
describes the non-magnetic mean field solution and weighted 
with 1 — Zcr. In the lower panel we have a similar comparison 
for [/ = 5 and /i = 0.15 concentrating on low energies. 
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FIG. 10: (Color online) A comparison of the imaginary parts 
of the transverse dynamic spin susceptibility for U — b and 
h — 0.0 (upper panel) and ^ = 0.15 (lower panel) calculated 
using renormalized perturbation theory (RPT, full line) and 
from a direct NRG calculation (dashed line). 



nary parts of the transverse susceptibility for h = 0.19, 
•we see that this overall agreement is maintained in the 
large field regime where we get metamagnetic behavior. 
The shapes of the low energy peaks for both quantities 
are well reproduced by the RPT formulae. Note that the 
peak in the real part is not at w = 0, so it is not fixed 
by the condition that determines t/t, but nevertheless is 
in good agreement with the NRG results. Due to their 
very small values it becomes difficult to calculate Za{h) 
as the system approaches localization for larger fields. In 
this regime as Za{h) the free quasiparticle density of 
states will converge to a delta-function. Self-energy cor- 
rections to the free quasiparticle propagators, which were 
used in the calculation of X(T.cr'{<-^), will become increas- 
ingly important as this Hmit is approached. Once the 
system has undergone the localization transition, and is 
completely polarized, however, we find that the values jia 



{za-{h) — 1) deduced from the self-energy give a quasipar- 
ticle density of states coinciding with the DMFT-NRG 
result of an upper and lower semi-circular bands. 

Results for the longitudinal susceptibility are shown in 
figures [12] and [131 In figure [12] we give the values for the 
real part as a function of uj for h = and h = 0.15. 
Here the peak height, which is at — 0, is fixed by 
the condition which determines C//. The widths of the 
peaks in the two sets of NRG results, however, are given 
reasonably well by the RPT equations. The imaginary 
part of the longitudinal susceptibility obtained by the two 
methods is given in figure [13] for h = 0.15. Again there 
is overall agreement between the two sets of results. The 
slight undulations seen in the RPT results are due to the 
sharp cut off in the band edges in the free quasiparticle 
density of states. For larger values of h the agreement 
with the NRG results is not as good as that as for the 



8 



60 



40 



3, 20^ 



-20 



'1 


' 1 ' 




Vfani /"DDT"! 




Imx (RPT) 




\ Rex (NRG) 




I Imx (NRG) _ 


: 








— 1 

— — VVi' 








,1 





-0.2 



-0.1 




CO 



0.1 



0.2 



FIG. 11: (Color online) The real and imaginary parts of the 
transverse dynamic spin susceptibility for U = 5 and h = 0.19. 
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FIG. 12: (Color online) The real part of the longitudinal dy- 
namic spin susceptibility for U = 5 and h = and h = 0.15. 
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FIG. 13: (Color online) The imaginary part of the longitudi- 
nal dynamic spin susceptibility for U = 5 and h = 0.15. 



transverse susceptibility, and the central peak in the real 
part of the RPT results narrows more rapidly with h than 
in those obtained from the direct NRG calculation. 



C. Strong coupling regime 

Finally we consider the strong coupling regime with 
U > Uc, where for h — the spectral density has a Mott- 
Hubbard gap so that for half-filling the system is an insu- 
lator. The electrons will be localized with free magnetic 
moments coupled by an effective antiferromagnetic ex- 
change J ^ /U . In fields such that h > J, the system 
polarizes completely as can be seen in figure fT4l where we 
show the total density of states p{uj) = p^{u}) + Pi{to) for 
h^Oandh = 0.2. 
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FIG. 14: (Color online) The total local spectral density p{uj) 
{oT U = 6 {ov h = (dashed line), Mott insulator, and h = 0.2 
(full line), fully polarized band insulator. 

For smaller fields, such that h < J, we do not find a 
convergent solution to the DMFT equations, and the it- 
erations oscillate between local states which are either 
completely full or empty. We interpret this as due to 
the tendency to antiferromagnetic order which in a weak 
field, due to the absence of anisotropy, will be almost 
perpendicular to the applied field in the x-y plane with a 
slight canting of the spins in the z-direction (spin fiopped 
phase). In this calculation no allowance has been made 
for this type of ordering, but this state can be well de- 
scribed using an effective Heisenberg model for the local- 
ized moments. 



IV. RESULTS AWAY FROM HALF FILLING 

After the extensive discussion of the behavior of the 
Hubbard model in a magnetic field at half filling, we 
want to compare these results with the situation where 
the system is doped with holes. As is well known, doping 
retains the metallic character of the system and one does 
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not find a paramagnetic metal-insulator transition any- 
more. Thus we do not find distinct regimes (a)-(c) any 
longer. To illustrate the characteristics of the magnetic 
response for the doped system we focus on two cases, one 
at quarter filling and one very close to half filling. 



A. Quarter Filled Case 

First we compare the results in the intermediate cou- 
pling regime with U = 5 at half-filling with those at 
quarter filling, x = 0.5. In the latter case the Fermi level 
falls in the lower Hubbard peak in the spectral density. 
To see how the band changes with increasing magnetic 
field we plot the density of states for both spin types, for 
the majority spin electrons in figure [15] (upper panel) and 
for the minority spin electrons in figure flSl (lower panel), 
for various values of the magnetic field. 
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FIG. 15: (Color online) The local spectral density for the 
majority spin Pt{(-o) (upper panel) and for the minority spin 
Pliuj) (lower panel) tor U — 5, x = 0.5 and various fields h. 



In the majority spin case the lower peak gains spectral 
weight on the low energy side and the weight in the up- 
per peaks decreases with increase of the field. For the 
fully polarized case {h > 0.4) the Fermi level, which is 



indicated by a dotted line, lies in the middle of the lower 
band, which has the non-interacting semi-circular shape. 
The opposite features can be seen in the minority spin 
case, with the spectral weight in the lower peak below 
the Fermi level decreasing and the weight in the upper 
peak increasing. Thus the increase of spectral weight be- 
low the Fermi level for the majority spin electrons, and 
the decrease for the minority spin electrons, can be seen 
to be due to a change of band shape rather than a sim- 
ple relative shift of the two bands, which would be the 
case in mean field theory. In the fully polarized state 
there are no minority states below the Fermi level and 
the upper peak in the majority state density of states 
has disappeared. 

The corresponding values for the inverse of the quasi- 
particle weight l/zcr{h) are shown in figure flGl for a range 
of fields. 




FIG. 16: (Color online) The inverse of the quasiparticle weight 
Zc{h) calculated from the impurity fixed point (FP) and di- 
rectly from the self-energy for U = 5, x = 0.5 and various 
fields h. The magnetization m{h) is also displayed. 



As noted in the impurity case^i, the quasiparticle weights 
differ for the two spin types with z^{h) > zi{h). The val- 
ues of Za-{h) have been calculated, as described earlier, 
both from the energy levels (FP) and from a numerical 
derivative NRG derived self-energy. There is reasonable 
agreement between the two sets of results, and the small 
differences can be attributed to the uncertainty in the nu- 
merical derivative of the NRG self-energy. As observed in 
the case for the non-symmetric Anderson model^i, there 
is an initial decrease of zi{h) with increase of h, whereas 
z^ih) increases monotonically. This implies that the ef- 
fective mass of the majority spin electrons decreases to 
its bare value, whilst the effective mass of the minor- 
ity spin electrons does not decrease much. The reason 
for that is that in the polarized system the up electrons 
cannot interact through the Hubbard interaction term, 
whereas a down spin electron can interact with all the 
up spin electrons leading to an enhanced mass. The field 
dependence of the magnetization is also shown in figure 
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[T6j and is similar to the half-filled case with a weak in- 
teraction {U = 2). We have calculated, but do not show, 
the corresponding occupation values for n° which again 
agree well with the values of fia-, confirming Luttinger's 
theorem in a magnetic field. 

We give two examples of results for the susceptibilities 
for this case. In figure flTl (upper panel) we plot the real 
and imaginary parts of the transverse susceptibility. 
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FIG. 17: (Color online) The real and imaginary parts of the 
transverse dynamic spin susceptibility (upper panel) and of 
the longitudinal dynamic spin susceptibility (lower panel) for 
[/ = 5, X = 0.5 and h = 0.1. 



Despite the large value of U, we can see that the peak 
heights are very much reduced compared with those seen 
in the half-filled case for U = 5. The peak widths are 
also an order of magnitude larger as can be seen from 
the w-scale. The RPT results reproduce well the overall 
features to be seen in the NRG results, but we note some 
discrepancies in the shape of the curve at larger frequen- 
cies, where the RPT shows more pronounced features. 
The real and imaginary parts for the longitudinal sus- 



ceptibility are shown in figure [17] (lower panel) . Again 
all the low energy features are reproduced in the RPT 
results and differences are mainly seen for tails at larger 
energies. In this regime, apart from the overall factor 
of 2, there is less difference between the transverse and 
longitudinal susceptibilities than at half-filling. 

Our conclusion from these results, and from calcula- 
tions with other values of U at quarter filling, is that 
when there is significant doping, the behavior in the 
field corresponds to a weakly correlated Fermi liquid, 
very similar to that at half-filling in the weak interaction 
regime. The only remarkable difference in the presence 
of a magnetic field is the spin dependence of the effective 
masses as shown in figure flGl 



B. Near half filling 

Very close to half-filling and for large values of U we 
have a qualitatively different parameter regime. Here the 
system is metaUic but we can expect strong correlation 
effects when U is of the order or greater than Uc, due 
to the much reduced phase space for quasiparticle scat- 
tering. We look at the case with 5% hole doping from 
half-fiUing and a value U = 6, which is just greater than 
the critical value for the metal-insulator transition. We 
show the spectral density of states for the majority spin 
state in figure [TSl (upper panel) and for the minority spin 
state in figure [H] (lower panel) for various values of the 
magnetic field. 

There is a clear sharp quasiparticle peak for h — aX the 
Fermi level (marked by a dotted line) at the top of the 
lower Hubbard band. As in the quarter fiUing case with 
U = 5, we see a similar transfer of spectral weight with 
increasing field to below the Fermi level for the major- 
ity spin case, and above the Fermi level for the minority 
spins. For large fields when the system is completely po- 
larized the Fermi level lies close to the top of the lower 
band in the majority spin spectrum. One can see in the 
lower panel that there is still a sharp narrow peak in 
the spectral density of the minority spin states above the 
Fermi level, though the spectrum for the majority states 
below the Fermi level is that of the non-interacting sys- 
tem. A spin up electron added above the Fermi level 
feels no interaction as the system is completely spin up 
polarized so these electrons see the non-interacting den- 
sity of states. On the other hand a spin down electron 
above the Fermi level interacts strongly with the sea of 
up spin electrons. The self-energy due to scattering with 
particle-hole pairs in the sea creates a distinct resonance 
in the down spin density of states just above the Fermi 
level. Just such a resonance was predicted by Hertz and 
Edwards^ for a Hubbard model in a strong ferromag- 
netic (fully polarized) state. 

The field dependence of the inverse of the quasiparticle 
weight is presented in figure [TH 

Again we find reasonable agreement between the two 
methods of calculation for these quantities. The mag- 
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FIG. 18: (Color online) The local spectral density for the 
majority spin Pi(tt') (upper panel) and for the minority spin 
Pl{uj) (lower panel) for U — G, x — 0.95 and various fields h. 
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FIG. 19: (Color online) The inverse of the quasiparticle weight 
Za-{h) calculated from the impurity fixed point (FP) and di- 
rectly from the self-energy for f/ — 6, x — 0.95 and various 
fields h. The inset shows the magnetization m{h). 



netization as a function of h is shown as an inset in the 
same figure. The behavior of Z|(/i) and zi{h) as a func- 
tion of h contrasts sharply with the behavior found for 
the metallic state at half-filling with U = b shown in fig- 
ure [6l For zero field the quasi-particle weight has a very 
similar value in both cases. At half-fiUing the tendency 
of the magnetic field to induce localization resulted in 
values of zi{h) and zi{h) {z^{h) = zi{h)) which decrease 
sharply as a function of h. In the 5% doped case with 
U — 6, the system must remain metallic and the quasi- 
particles weights, z^{h) and zi{h), both increase in large 
fields though their values differ significantly. The quasi- 
particle weight for the minority spin electrons decreases 
initially with increase of h, whereas that for the majority 
spins zi{h) increases monotonically and quite dramati- 
cally with h. When the system becomes fully polarized 
{h ~ 0.26) the up spin electrons become essentially non- 
interacting, z-i{h) = 1, whereas there is a strong inter- 
action for a down spin electron and we find in this case 
zi{h) ~ 0.15. The interpretation for this is as given in the 
previous paragraph for the spectral densities. On further 
increasing the magnetic field zi{h) also tends to one, but 
relatively slowly, as can be seen in figure [T9l Note that 
the results in this regime are based on the calculation 
from the self-energy, as the method based on the fixed 
point analysis becomes difficult to apply in this regime. 

Laloux et ali^ compared the quasiparticle weight at 
half filling for the infinite dimensional model with results 
from the Gutzwiller approximation. The values for the 
infinite dimensional model were found to be significantly 
smaller than the Gutzwiller predictions, the ratio is more 
than a factor of 2 for [/ > 4 and zero field. Spalek et 

4,28 iiave made predictions based on the Gutzwiller ap- 
proach for situations away from half filling in finite field, 
Z| zi. As in the study of Laloux et alr^ our results for 
Zi and zi are significantly smaller than the Gutzwiller 
predictions. 

For the fully polarized case {h = 0.26) we show the 
comparison of the weighted free quasiparticle density of 
states ZcPo,cr{^) with the full spectrum Pai^) in figure 

m 

Note that the parameters flo,a- and z^ in pQ^a-{uj) are 
purely derived from the NRG self-energy in this case. 
We can see that the different values for the field depen- 
dent quasiparticle weight Za-{h) for up and down spin lead 
to remarkably different quasiparticle band shapes. With 
Z| ~ 1 the majority spin quasiparticle band is essentially 
that of the non-interacting density of states. The very 
much smaller value Z| leads to a narrow quasiparticle 
band above the Fermi level. The low energy fiank of this 
quasiparticle band describes well the narrow peak seen 
in the spectral density just above the Fermi level. To de- 
scribe these strong asymmetries in the spectral densities 
near half filling, we need Z| ^ ^x, which contrasts with 
the cases at half filling such as in figure [9] where always 
= z^. This suggests a discontinuous behavior of the 
renormalization factors Za on the approach to half fiUing. 
Also for this case we show plots for the two suscepti- 



12 




10 



FIG. 20: (Color online) The free quasiparticle density of states 
in comparison with interacting local spectral density for U = 
6, s = 0.95 and h = 0.26. 



bilities for a field oi h — 0.15. In figure [2T] (upper panel) 
we give the real and imaginary parts of the transverse 
susceptibility. 

The low energy features are seen on an a;-scale an order 
of magnitude smaller than that for quarter filling due to 
the much stronger renormalization effects in this regime. 
There is excellent agreement both with the peak positions 
and shapes between the NRG and RPT results for both 
quantities. This is also seen to be the case for the real and 
imaginary part of the longitudinal susceptibility shown in 
figure ED (lower panel) , though the peak in the real part 
can be seen to be marginally narrower in the RPT results. 

At the end of this section, we conclude that already a 
small doping of the system is enough to maintain a metal- 
lic character even for very strong interaction. Although 
the zero field spectra of the half filled case for U — 5 and 
the small doping case with U = 6 display very similar 
zero field behavior, i.e. a strongly renormalized quasipar- 
ticle band with similar Za, no field induced localization 
transition occurs for finite doping and no metamagnetic 
behavior is observed in the latter case. 



V. SUMMARY 

In this paper we have used the DMFT-NRG method to 
calculate the spectral densities for one-particle and two- 
particle response functions for the infinite dimensional 
Hubbard model in a magnetic field, for the qualitatively 
different filling regimes and interaction strengths. The 
results extend earlier calculations of Laloux et al.- using 
the ED method, which were restricted to the case of half- 
filling. Our results there are on the whole consistent with 
this earlier work, except in the insulating regime for weak 
fields, where we could not find a convergent solution of 
the DMFT equations. We attributed this to the fact 
that in this regime the magnetic field is smaller than the 
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FIG. 21: (Color online) The real and imaginary parts of the 
transverse dynamic spin susceptibility (upper panel) and of 
the longitudinal dynamic spin susceptibility (lower panel) for 
U = 6,x^ 0.95 and h = 0.15. 



exchange coupling between the localized spins so that 
the ground state would be one in which the spins would 
have a canted antiferromagnetic ordering in the plane 
perpendicular to the field. 

Well away from half filling we find a magnetic re- 
sponse similar to the weakly correlated case even for 
large values of U. The large phase space for quasipar- 
ticle scattering in this regime leads to modest renormal- 
ization effects. Here, we find spin dependent quasiparti- 
cle weights, z^{h) ^ zi{h). This implies spin dependent 
as well as field dependent effective masses, which have 
been discussed earlier in work by Spalek et and 
Riseborough~. The calculations by Spalek et al. were 
based on a Gutzwiller~ and a mean-field slave boson 
approach^. We can make a comparison of our results 
(section VI. B) near half filling, x — 0.95, with theirs in 
the later work<4. We find a qualitatively similar behavior 
with the majority spin effective mass decreasing with /i, 
but quantitatively there are differences. The field depen- 
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dence of the minority spin effective mass shows 
a very slow increase initially in both sets of results, but 
the large field behavior is quite different. As seen in fig- 
ure [19] we find a significant decrease in for large 
fields whereas the corresponding quantity in figure 3 in 
reference! increases. 

The strong magnetic field dependence of the effective 
masses found in the calculations by Riseborough is based 
on the assumption that the system is close to a ferromag- 
netic transition (paramagnon theory). However, DMFT 
calculations for the Hubbard model find that any ferro- 
magnetism in the Hubbard model only occurs in a very 
small region of the parameter space near half-filling and 
for very large values of /7.— Our results are well away 
from this regime and the large effective masses obtained 
here can be attributed to the tendency to localization 
rather than the tendency to ferromagnetism. 

Using the field dependent renormaHzed parameters 
Zcr{h) and jlo,^{h) in the RPT formulae for the dynamic 
local longitudinal and transverse spin susceptibilities we 
found agreement with the overall features to be seen in 
the DMFT-NRG results for these quantities. In the case 
of the transverse spin susceptibility excellent agreement 
was found in all the metallic regimes and for all val- 
ues of the magnetic field considered, except in the high 
field regime at half filling as the localization point is ap- 
proached, where consistent values of the renormalized 
parameters are difficult to calculate. The comparison 
of the RPT results with those from NRG was also ex- 
cellent for the longitudinal dynamic susceptibility in the 
weaker field regime h < 0.15 but less good for higher 
fields, h > 0.15. 

In all metallic parameter regimes a spin dependent 
Luttinger theorem in the form Ua- = n'^, the number of 
particles equals the number of quasiparticles, was found 
to be satisfied for all strengths of the magnetic field. In 
this form it even holds in the fully polarized insulating 
state. 



Phenomena like field and spin dependent effec- 
tive masses and metamagnetic behavior have been 
observed experimentally in several heavy fermion 
compounds^i^i^i^. The Hubbard model, however, being 
a one band model is not an appropriate starting point to 
make a quantitative comparison with the heavy fermion 
class of materials. A periodic Anderson model with a two 
band structure and including the degeneracy of the f elec- 
trons would be a better model to describe these materials. 
Field dependent effects in this model have been studied 
by several techniques, modified perturbation theory^^, 
exact diagonalization2S,, 1/N expansion^ and variational 
approach^. The approach used here could be general- 
ized to the periodic Anderson model, but restricted to 
the non-degenerate case and TV = 2 as it is computa- 
tionally too demanding in the NRG to deal with higher 
degeneracy. The Hubbard model at half filling has been 
used as a lattice model to describe the strongly renor- 
malized Fermi liquid '^H&&2. However, the metamagnetic 
behavior predicted for relevant parameter regime is not 
seen experimentallySfi. In section VI. B we found for small 
doping large effective masses, but no metamagnetic be- 
havior. This raises the possibility that the weakly doped 
Hubbard model could serve as a basis for interpreting the 
experimental results for Hquid '^He. 



Acknowledgment 

We wish to thank N. Dupuis, D.M. Edwards, W. Roller, 
D. Meyer and A. Oguri for helpful discussions and W. 
Roller and D. Meyer for their contributions to the devel- 
opment of the NRG programs. One of us (J.B.) thanks 
the Gottlieb Daimler and Rarl Benz Foundation, the Ger- 
man Academic Exchange Service (DAAD) and the EP- 
SRC for financial support. 



^ W. Joss, J. M. van Ruitenbeek, G. W. Crabtree, J. L. 

Tholence, A. P. J. van Deursen, and Z. Fisk, Phys. Rev. 

Lett. 59, 1609 (1987). 
^ R. G. Goodrich, N. Harrison, A. Teklu, D. Young, and 

Z. Fisk, Phys. Rev. Lett. 82, 3669 (1999). 
^ H. Aoki, S. Uji, A. K. Albessard, and Y. Onuki, Phys. Rev. 

Lett. 71, 2110 (1993). 

P. Korbel, J. Spalek, W. Wojcik, and M. Acquarone, Phys. 

Rev. B 52, R2213 (1995). 
^ M. Manekar, S. Chaudhary, M. K. Chattopadhyay, K. J. 

Singh, S. B. Roy, and P. Chaddah, J. Phys.: Cond. Mat. 

12, 9645 (2000). 
^ D. Vollhardt, Rev. Mod. Phys. 56, 99 (1984). 

L. Laloux, A. Georges, and W. Krauth, Phys. Rev. B 50, 

3092 (1994). 

* V. Janis and G. Czycholl, Phys. Rev. B 61, 9875 (2000). 
® F. Kagawa, T. Itou, K. Miyagawa, and K. Kanoda, Phys. 
Rev. Lett. 93, 127001 (2004). 



" A. C. Hewson, J. Bauer, and W. Roller, Phys. Rev. B 73, 
045117 (2006). 

" J. Bauer and A. C. Hewson (2007), cond-mat/0705.3818. 
12 A. C. Hewson, Phys. Rev. Lett. 70, 4007 (1993). 
" A. C. Hewson, J. Phys.: Cond. Mat. 13, 10011 (2001). 
" A. C. Hewson, J. Phys.: Cond. Mat. 18, 1815 (2006). 
1^ W. Metzner and D. Vollhardt, Phys. Rev. Lett. 62, 324 
(1989). 

E. Miiller-Hartmann, Z. Phys. B 74, 507 (1989). 

A. Georges, G. Kotliar, W. Krauth, and M. Rozenberg, 

Rev. Mod. Phys. 68, 13 (1996). 
1** J. M. Luttinger, Phys. Rev. 119, 1153 (1960). 
" H. R. Krishna-murthy, J. W. Wilkins, and K. G. Wilson, 

Phys. Rev. B 21, 1003 (1980). 
2° A. C. Hewson, A. Oguri, and D. Meyer, Eur. Phys. J. B 

40, 177 (2004). 

R. Bulla, A. C. Hewson, and T. Pruschke, J. Phys.: Cond. 
Mat. 10, 8365 (1998). 



14 



R. Peters, T. Pruschke, and F. B. Anders, Phys. Rev. B 
74, 245114 (2006). 

A. Weichselbaum and J. von Delft (2006), cond- 
mat/0607497. 

F. B. Anders and A. Schiller, Phys. Rev. Lett. 95, 196801 
(2005). 

R. Bulla, Phys. Rev. Lett. 83, 136 (1999). 

J. Bauer, A. C. Hewson, and A. Oguri, J. Magn. Magn. 

Mat. 310, 1133 (2007). 

J. A. Hertz and D. M. Edwards, Phys. Rev. Lett. 28, 1334 
(1972). 

J. Spalek and P. Gopalan, Phys. Rev. Lett. 64, 2823 
(1990). 

P. S. Riseborough, Phil. Mag. 86, 2581 (2006). 



R. Zitzler, T. Pruschke, and R. Bulla, Eur. Phys. J. B 27, 
473 (2002). 

S. V. Dordevic, K. S. D. Beach, N. Takeda, Y. .J. Wang, 
M. B. Maple, and D. N. Basov, Phys. Rev. Lett. 96, 017403 
(2006). 

D. Meyer and W. Nolting, Phys. Rev. B 64, 052402 (2001). 

T. Saso and M. Itoh, Phys. Rev. B 53, 6877 (1996). 

Y. Ono, ,J. Phys. Soc. Japan 67, 2197 (1998). 

D. Edwards and A. C. M. Green, Z. Phys. B 103, 243 

(1997). 

O. Buu, A. Forbes, A. van Steenbergen, S. Wiegers, 
G. Remtoyi, L. Puech, and P. Wolf, J. Low Temp. Phys. 
110, 311 (1998). 



